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We have applied the Andreev-Lifshitz hydrodynamic theory of supersolids to an ordinary solid. 
This theory includes an internal pressure P, distinct from the applied pressure P a and the stress 
tensor Xik- Under uniform static P a , we have Xik = (P~ P a )Sik- For P a ^ 0, Maxwell relations imply 
that P ~ P^ . The theory also permits vacancy diffusion but treats vacancies as conserved. It gives 
three sets of propagating elastic modes; it also gives two diffusive modes, one largely of entropy 
density and one largely of vacancy density (or, more generally, defect density). For the vacancy 
diffusion mode (or, equivalently, the lattice diffusion mode) the vacancies behave like a fluid within 
the solid, with the deviations of internal pressure associated with density changes nearly canceling 
the deviations of stress associated with strain. We briefly consider pressurization experiments in 
solid 4 He at low temperatures in light of this lattice diffusion mode, which for small P a has diffusion 
constant Dl ~ Pa- The general principles of the theory - that both volume and strain should be 
included as thermodynamic variables, with the result that both P and Xik appear - should apply 
to all solids under pressure, especially near the solid-liquid transition. The lattice diffusion mode 
provides an additional degree of freedom that may permit surfaces with different surface treatments 
to generate different responses in the bulk. 

PACS numbers: 67.80. B-, 67.80.bd, 05.70.Ln, 63.10.+a 



I. INTRODUCTION 

Since the late 1960's there have been theoretical sug- 
gestions that solids might display flow behavior simi- 
lar to what is found in supernuidsj 1 - - — For that reason 
there has been a great deal of interest in solid 4 He as a 
candidate super solidM* The first experimental indication 
of superflow was the appearance of a non-classical mo- 
ment of inertia (NCRI), first observed by Chan's group, 
since confirmed by many other laboratories, and strongly 
linked to disorder^— In addition, the shear modulus 
shows anomalous behavior, 16 although not enough to ex- 
plain the NCRI experiments^ 7 - Non-NCRI superflow has 
been searched for but not observed We also note re- 
cent experiments that argue against any supersolid signa- 
ture above approximately 55 mK. 19 Further works cast- 
ing doubt on supersolidity are a study of bcc 4 He that 
shows unusual NCRI behavior at higher temperatures^ 
and a study showing that the NCRI behavior due to plas- 
ticity has different properties than due to quenching^ 1 - 

In a recent experiment on a pancake-shaped sample, a 
capacitance gauge monitored the pressure as a function 
of temperature Ti22 Samples were produced by both the 
slow-cooling blocked capillary method and by the more 
rapid quench-cooling method, which gives more disor- 
dered samples. In one set of measurements the sample 
was quench-cooled below 1 K in 144 s, during which 
time the pressure decreased. This is perhaps an indi- 
cation that vacancies, formed during the quench, were 
leaving the sample. For a blocked-capillary sample the 
temperature was lowered below 500 mK while the pres- 
sure was monitored vs T. The sample was then annealed 
at 1.65 K, where the pressure increased, perhaps an in- 
dication that vacancies now were entering the sample. A 
second cooldown yielded, by a reduced T 2 term in the 



pressure, an indication that the sample was less disor- 
dered, but that disorder remained. Even at a constant 
temperature of 19 mK the pressure continued to relax, 
which is consistent with vacancies equilibrating. The fact 
that the observed relaxation times do not saturate at the 
temperatures studied indicates that the temperature is 
not yet low enough that quantum relaxation processes 
dominate thermal relaxation processes. 

This experiment can perhaps be interpreted under the 
assumption that the system is not supersolid. We have 
therefore undertaken a theoretical study of the macro- 
scopic flow properties of a one-component ordinary solid. 
Our basis is the theory of Andreev and Lifshitz (AL) 
for the macroscopic behavior of a supersolid. They in- 
cluded volume V as an extensive variable, in addition to 
Wik = Vwik, where Wik is the non-symmetrized strain. 
This permitted them to continuously go to the super- 
fluid limit as Wik becomes irrelevant. The point of the 
present work is that on eliminating the superfluid vari- 
ables, the theory should apply to an ordinary solid. 1 We 
employ a variation on the notation of Ref. which gives 
a more explicit derivation of the equations of motion and 
extends Ref. [l| to include nonlinear termsj 24 ' 25 Note also 
the theory of Fleming and Cohen 2 ^ for an ordinary solid, 
which gives equations with a similar structure, and sim- 
ilar modes, but uses a very different notation (and does 
not consider an applied pressure P a ). Both Ref. |l| and 
Ref. |26| implicitly assume that uniform vacancy number- 
changing bulk processes are negligible, and neglect in- 
terstitials and impurities. Recently Yoo and Dorsey 27 
considered the effect of a lattice diffusion mode on light 
scattering by a supersolid, but also briefly considering an 
ordinary solid. 

As noted by Martin, Parodi, and Pershan^ the nor- 
mal system has eight degrees of freedom, given by two 
scalar thermodynamic quantities (which can be taken to 
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be the mass density p and the entropy density s) and two 
vector quantities: the lattice vector Ui and the velocity x>i 
associated with the momentum density gi = pvi. (With 
1714 the atomic mass and n the number density of 4 He 
atoms, we have p = m^n.) As a consequence there are 
eight normal modes. For a uniform infinite system these 
modes are three pairs of propagating elastic waves and 
two diffusive modes, one primarily of the temperature T 
and the other primarily of SjUj. In the absence of lattice 
defects, for a variation 5ui the relationship 

di(6ui) » -Sp/p (1) 

holds, giving the system one fewer degree of freedom, 
and thus one fewer mode. One can think of this missing 
mode, associated with the dynamical violation of JTJ, as 
being associated with vacancies, as noted in Ref. l28l 

The present work obtains the diffusion constant and 
the physical properties of this diffusive mode, for both 
zero and non-zero P a . ( 4 He must be under P a « 25 at- 
mospheres to solidify.) We find that the physical charac- 
ter of the mode is that it involves essentially zero stress 
deviation, because the fluid-like stress (associated with 
changes in mass density) nearly cancels the solid-like 
stress (associated with changes in strain). 

Allowing vacancies to move permits mass change with- 
out lattice motion^ This allows one to take the fluid 
limit of zero crystallinity, and study the evolution of the 
sound velocity as the system evolves from the perfect 
solid to perfect liquid. By perfect solid we mean one with 
no defects and a one-to-one relationship between lattice 
points and atoms; by perfect liquid we mean one with no 
lattice structure or, equivalently, one with no sensitivity 
to an imaginary lattice structure. A gel has properties of 
both, but is multi-component £L 

Section HD gives the form of AL supersolid theory when 
restricted to a normal solid, including the possibility of 
lattice defects. Although we specifically have vacancies in 
mind, 3 He impurities could be accounted for if its density 
were included as an additional thermodynamic variable, 
which would require extension of the AL theory. Note 
also the case of (two-component) superionic conductors, 
which includes certain high-temperature alkali halides, 
where the larger halide ions remain in a lattice but the 
lattice of the smaller alkali ions "melts." Section [TTT] dis- 
cusses elasticity and internal pressure for a crystal under 
static and uniform applied pressure P a , and calculates 
internal pressure P and strain. We find that P ~ P 2 , so 
that for small P a the effect of P is very small; see eq. (|34p . 
For small P a the strain is largely linear in P a , as expected, 
but there is a P 2 correction. Section IIVI derives the nor- 
mal modes for the ordinary solid. Section [V] considers 
how such modes can be generated (including the possi- 
ble effect of different surface treatments), and applies the 
theory to the pressurization experiments. 22 Section IVII 
provides a summary and our conclusions. Appendix [X] 
gives the thermodynamics and dynamics of the AL the- 
ory for the supersolid. Appendix IB"! calculates some ther- 
modynamic derivatives that appear in the normal modes 



in terms of P a . 



II. ANDREEV-LIFSHITZ NORMAL SOLID 
WITH DEFECTS 

We employ the primary quantities energy density e, 
lattice displacement Ui, and non-symmetrized strain 

w ik = diU k . (2) 

We consider a normal solid by setting p s = 0, p n = p. 
v n = v, and eliminating the superfluid equation from the 
equations for the supersolid (given in Appendix [A| . 

A. Thermodynamics 



The appropriate thermodynamic equations are 

de = Tds + Xikdwik + pdp + v ■ dg, (3) 

e = -P + Ts + \ lk Wik + pp + v- g, (4) 

= —dP + sdT + WikdXik + pdp + g ■ dv. (5) 



Here Xik of AL is an elastic tensor density (with units of 
pressure P), and p. is the chemical potential (with units 
of velocity squared); Xik is the same as ctq. of Ref. [3ll 

B. Dynamics 

The appropriate linearized equations of motion for the 
independent variables s, Ui, p, and t>, are 

d t s + difi = 0, (6) 

dm = Ui, (7) 

dtp + di9i = 0, (8) 

d t9l + d k n tk = 0, (9) 

where the fluxes fi (of entropy), (of momentum), gi 
(of mass), and the "source" Ui (terminology introduced 
here) are given by 

fi = sv i -'^d j T-^d l X lj , (10) 

Ui = v i + -2-d j T + ij d l \ lj , (II) 

ILfc = (PSih - X ik ) - mUmdmVl, (12) 

9i = pvi- (13) 

AL use Gik ~ —Hik^ The term in (TTT1) proportional to 
(3ij allows the lattice velocity £tj to differ from the velocity 
Vi associated with mass flow. It leads, as we show, to a 
lattice diffusion mode for which iij ^ Vi and neither is 
zero. 

Both the otij and terms can be rewritten as flux 
terms. Linearizing about equilibrium, with primes de- 
noting deviations from equilibrium, yields 

d t Ui + OjSij = v'i, (14) 
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where 

Si^-^T'-PuXjt, (15) 

In (fT4|) , Vi can be thought of as a "lattice source" , and 
Sij as a "lattice flux." The ftu term gives, in principle, 
anisotropic vacancy diffusion. 

Recall that a diffusion constant D is proportional to 
a characteristic velocity times a characteristic mean-free 
path, so it has units of m 2 /sec. In terms of a D, the 
dissipative coefficients have the following units: Kij has 
units of s times D; aij has units of D\ /3y has units of 
inverse pressure times D; and r\ikim has units of p times 
D. 



As employed by Ref. U this pressure term, in contrast 
to Xik alone (Ref. HH does not include P), permits one 
to continuously approach the superfluid limit, when the 
lattice disappears. In the present case, it permits one to 
continuously approach the ordinary liquid limit. 

The consequences of a nonzero P include, but are not 
limited to, a mode where vacancies are permitted to dif- 
fuse. Thermodynamic derivatives of P are essential for 
defect diffusion, and also affect the elastic modes. More- 
over, they are needed to obtain the pure liquid limit for 
longitudinal sound on letting the crystallinity go to zero. 
We first use a Maxwell relation to find an explicit expres- 
sion for P as a function of strain. 



III. CRYSTAL UNDER PRESSURE 

A. Internal Pressure and Elasticity 

The momentum conservation equation (0) implicitly 
contains the term Xik — P&ik , which determines the force 
on the surface of the solid. An internal pressure P does 
not appear in the thermodynamics of Ref. 3l|, which 
does not consider either a lattice under applied pressure 
or the presence of defects. However, the extensive en- 
ergy E = eV, which depends on the extensive variables 
(S, V, N, Wik = Vwik,Vg), has second derivatives that 
satisfy the Maxwell relation 



B. Internal Pressure P 



dP 



dXjk 

dW lk ~ dV ' 



(16) 



where the appropriate variables are held constant. For 
solid 4 He under an applied pressure P a , this makes P 
non-zero. 

In principle we may let E depend on the number 
of vacancies Ny, with associated "chemical potential" 
4>v = dE/dNy (with units of energy, rather than veloc- 
ity squared). Then the additional Maxwell relation 



dP 



dN v dV 



(17) 



follows, with the appropriate variables held constant. If 
the vacancies are not in equilibrium (i.e., (f>v 0), this 
also makes P non-zero. The general results of the present 
work (e.g., a nonzero lattice diffusion constant) can thus 
be made applicable to a solid not under P a but having 
vacancies out of local thermal equilibrium. Terms found 
here to depend on P a may in that case depend on the dif- 
ference between actual concentration of vacancies and the 
equilibrium concentration of vacancies. However, we ex- 
pect a P a of ~ 25 atm to dominate the effect of vacancies, 
and thus we neglect their effect on P. Although Refs. [0, 
[H and [H introduce the internal pressure P, they do 
not calculate P or its thermodynamic derivatives. Ref. U 
and the present work neglect the possibility of interstitial 
atoms^ For a reference that considers interstitials, see 
Ref. HI 



Since holding (V, N) constant is equivalent to hold- 
ing (V, p) constant, and similarly for (S, N) and (<r = 
s/p,N), we use these sets interchangeably. We rewrite 
© as 



dP 



dW lk 



1 dP 



dXik 



V,(T, P 



dV 



V,S.N V dw ik 

For constant Wik we have 

= dWik = w lk dV + Vdwik, 



W ik ,S,N 



so that 



Then 

dXik 



dwi 



ik 



dV 



W ik ,S,N 



dX t 



f)V 

uv W zkl S,N 

and (fT8l) gives 
dP 



dV 



Wjk 

' V ' 



Wji dX lk 



(18) 



(19) 



(20) 



(21) 



V,a,p 



dwik 



= -v 



dXik 



V,a,p 



dV 



= _y 9X V 



dV 



W ik ,S,N 



dXik 



dw 



'jl 



(22) 



V,CT, P 



We employ Ref. [3l| for the elasticity tensor Xik in an 
isotropic solid. Using superscript (0) to denote the equi- 
librium value of Xik and the strain u;^, we have 

Ak = (k - ^4 0) + ^ (4 0) + »2 J ) ■ (23) 

where K and p,y are the bulk and shear moduli, and both 
and uk? are to be determined under a given applied 
pressure P a . Eq. (|22|) then gives 

\ w ik + w ki ) i 

(24) 



dP 



dw ik 



V,CT, P 



= \K* - ) Sikwffl 
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where 



K* =K -V 



a result that appears to be new. Further, A^ = S \k , 
where 



.(0) 



dK 
dV 



Wiu .<7.N 



dV 



Wik ,<T,N 

(25) 



Under uniform P a we expect an isotropic response, so 
(0) 



Then ([51} becomes 
dP 



dwik 



Oik (0) 

3 " 



(26) 



(27) 



V,a,p 



Integration of (|2~T1) with respect to Wik gives the part of 
the internal pressure dependent on the strain to be 



P = 



(28) 



where we take the integration constant to be zero^i For 



w ii = w 22 = w 33 > we then have 



(0)' 



P = -K*w' n 



(29) 



This result applies to the case of a strongly crystalline 
material. In the opposite limit where the crystallinity 
disappears and the particles are weakly interacting, part 
of P would be given by the ideal gas law. 



C. Strain Wik 

As discussed above, under an applied pressure the force 
on the surface of a solid is 



ik • 



Taking the trace yields 



A 



CO) 
ll 

3 



P = ~Pn 



Substitution from ([23} and ([29} gives 



3Kw^ - \K*wfy = -P a . 



(30) 



(31) 



(32) 



Since an applied pressure should cause a negative strain, 
only the solution for w^i < is physical. 

For solid 4 He, we expect both and P a /K to be 
small. The solution of ([321) to second order in P a is 



(33) 



,.(0) ~ P « ■ P a K * 
11 ~ 3K 6K 3 



The first term is what one would get on neglecting P in 
([3T} . To second order in P a /K, eq. (|29f) then gives 



A£? = -P a + K* Pa 



2K 2 ' 



(35) 



The first term in Ajj is what one obtains on neglecting 
P in ([5T} , and in agreement with Ref . [3l|. 



IV. NORMAL MODES OF ANDREEV-LIFSHITZ 
NORMAL SOLID WITH DEFECTS 



As noted earlier, this system has eight variables: s, p, 
gi and it j . Disturbances from equilibrium will be denoted 
by primes, so we use s', p' ', g[ « pv[, and m^. There are 
correspondingly eight normal modes. For an infinite sys- 
tem we assume a disturbance of the form exp[i(k-f—uit)], 
where the real wavevector k is considered to be known, 
but U) is unknown. For the disturbance to decay in time, 
Im(uj) < 0. Six modes come in three degenerate pairs, 
with g[ and u\ strongly coupled, and correspond to ordi- 
nary elasticity. The other two modes are diffusive, with 
temperature diffusion nearly decoupled from lattice diffu- 
sion. To ensure this decoupling we set the (off-diagonal) 
temperature-lattice transport coefficient onj — 0, and set 
the distinct but similar-looking thermal expansion coeffi- 
cient a = i 31 ' 36 We consider an isotropic solid, for which 
Kij = nSij and /3y = /3#y , this /? not to be confused with 
the identical symbol sometimes used for the thermal ex- 
pansion coefficient. 37 

We also neglect the tensor viscosity r)iki m > which to 
lowest order in k does not contribute to the modes. The 
fluctuation of the tensor Il- fe ([T2"} has a term from the 
viscosity ~ r)ikimk m Vi and a term from the stress tensor 



A; 



i k 



(dXik/dw^w' Then, by ©, A- 



ik 



kjv'i/u. Thus, for both propagating modes (u ~ k) and 
diffusive modes (ui ~ fc 2 ), the term in Ii' ik due to viscosity 
is, at the least, of order k relative to the term X' ik , and is 
therefore neglected in the long wavelength limit. 



A. Thermal Diffusion 

For the normal solid it is convenient to work with p 
and a — s/p as variables, because a diffuses but does 
not flow, and therefore is nearly conserved. To see this 
note that, to lowest order in deviations from equilibrium, 
eq. © and flgj) yield 



T(da/dT) p 



W, (36) 



where we have used a — 0. 3S This equation describes 
entropy diffusion, with a' ^ and 



P = K* 



2K 2 ' 



(34) 



oj = —iDxk , Dq 



P T(da/dT)^ 



(37) 
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For this mode u\ = v[ — p' = 0. If a is small but non- 
zero the frequency will not change to lowest order in a, 
but from the equations for p, g, and u these quantities 
would develop amplitudes proportional to a' and a, and 
thus have negligible amplitude as a — > 0. We consider 
only the case where the effects of a can be neglected. 



B. Elastic Modes 

We obtain the elastic modes by taking a' = and ne- 
glecting dissipative and nonlinear terms in ([T])-©. Thus, 
eq. |7]) gives u[ = v[. In the remainder of this work, all 
thermodynamic derivatives are taken at constant a, and 
derivatives with respect to p are taken at constant Wi k 
and vice-versa, unless otherwise specified. Further, when 
derivatives with respect to a specific component of Wi k 
are taken, the other components of Wi k are held fixed. 
Then by $TZ§ and (13]), eqs. © and © become^ 



dP 

= pu'i + -g^dip' + 
= p' + pdiu'i. 



dP 

Own 



dr 



dp 



dkp' 



Own 



dkw'ji, 
(38) 
(39) 



Clearly, a' does not couple to the other variables. On 
linearizing, eq. (|39|) gives p' = —pdiu'^ so with (f2|), 
eq. ([38]) becomes 



dP 

= pu'i - p—did k u' k 



dP 



dwji 
dX ik i 
p ~dp~ dkd i u i 



dX 



dw 



•ji 



d k d 3 u[. (40) 



The second term gives the pure fluidlike (longitudinal) 
response, which occurs for P ^ (e.g., an imperfect solid 
or a solid under P a ), and the fifth term gives the pure 
solidlike (longitudinal and transverse) response. 

Appendix [Bl shows that, for uniform static P a , certain 
quantities are isotropic. This permits us to define 



dP dP 



dw 



3l 



dw ' 



dXji = dX _ ^ 

dp dp 3 dw 



(41) 



Appendix [Bl also shows that 
dX ik 



dw 



ji 



dX 

g-SikSji + p v (StjSki + 5 k j5a - 2SikSji) . 



(42) 



Thus (gOD gives 



dP 

pu'i - p-^did k u' k 



P ~dp 1 kUk 



OX 
dw 



dP 

-^-did k u' k 
dw 

- pv ) didju'j 



MvV 2 <. (43) 



On letting di iki and dt —> —icu, eq. (|43|) becomes 



fa (-puj 2 + pv)^^ 
dP dP 



dp di 



dX 
dp 



dX 

dw 



Pv 



h{k-u). (44) 



Longitudinal Mode: If k ■ u ^ 0, then (j44|) shows 
that Ui is along fej, so the mode is longitudinal. Moreover, 
eq. (|44[) gives the normal mode frequencies 



u 2 = 



dP 
dp 



1 dP 

p dw 

-is] k 2 



dx 

dp 
cfk 



1 dX 
p dw 



(45) 



where 



HL 



dP dX 2 
Jp~ _ d~p' ° ls 



IdX IdP 

-o o-- 46 

p dw p dw 



The liquid-like velocity c;l contains thermodynamic 
derivatives with respect to the density p, and the solid- 
like velocity cis contains thermodynamic derivatives with 
respect to the strain Wik ■ Eq. (|45|) gives a velocity for lon- 
gitudinal sound that is similar to that found in Ref. 128 . 

Appendix IB1 finds the four derivatives in (|46[) in terms 
of P a , which to second order in P a /K give 



Pn (K 



1L 



C IS 



P 2 K* ( K* p dK* 
1 - 



p \ K 1 J ' 2pK 2 V~ K ' K* dp 

K+lp v PaK*_ _ P 2 K* 2 
p p K 2pK 3 



(47) 
(48) 



where K* is defined in (|25p. For P a — we have cf = 
cf L + c 2 s = [K + (4/3)/zy]/p, which agrees with Ref. HH 
for an ordinary solid. 

Transverse Mode: If k ■ u = 0, so that the mode is 
transverse, then (|44[) gives the normal mode frequencies 



j 2 = ^k 2 . 



(49) 



P 

From (f3"9"]l . for the transverse mode /?' = 0. Eq 
agrees with Ref. [3l] for an ordinary solid. 

For both longitudinal and transverse mode frequencies 
eq. (|3l))) is satisfied by a' = 0. 



C. Lattice Diffusion 

The lattice diffusion mode is the most subtle of the 
modes. For this mode, as for the elastic modes, we con- 
sider that a is constant, but we do not take v[ = u^. 
Rather, we assume that to = —iDi,k 2 , where the lattice 
mode diffusion constant > is to be determined, and 
we keep the dissipative terms in the equations of motion 
for v[, u'i, and p' . 
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With = (35ij (i.e., an isotropic solid), onj 
setting a" = 0, eqs. (JSJ and (J7J give 



-tup 



v' i + p(ikk)Kk- 



0, and 

(50) 
(51) 



If we assume that the mode is longitudinal, with v[ ~ fcj 
(the consistency of this assumption to be determined 
below), then the first of these equations implies that 
v'i ~ kip'. Therefore in ((9|) the term dtgi ~ wpv[ is of 
order k 2 relative to the hp' dependence of dk^l' ik , and 
is neglected in the long wavelength limit. As a conse- 
quence, d\jV ik ~ 0: the contributions from the liquid-like 
part P'Sik and from the solid-like part — X' ik nearly can- 
cel. This can only occur for an imperfect solid or a solid 
under applied pressure P a , which has both liquid- like and 
solid-like responses. 

Thus, neglecting the d t gi ~ uipv^ term and neglecting 
the viscosity r\iUm (as discussed above), eq. © gives 



ikip — 
dp 



, dX lk dX 



dp dwji 
Substitution from flU) and gives 



ik I 
U, 



dP 



dw 



(52) 



,dP 
dp 



ikip 



,dX 
dp 



dX_ 
dw V 



dP 

kihu', - pvk 2 u'i + —kikiu'i. (53) 
dw 



All but one term in (|53l) is along fe,, and the remaining 
term is along u[. Therefore we deduce that itj is along 
hi, and thus kikiu', = fc 2 w^. Then (|53| becomes 



ip 



— - —\ - -k ' (— - — 

dp dp J 1 1 \dw dw 



(54) 



Further, eq. 

(dXik/dwjfiw'^ - 

ki , 



(fSTj) gives, on taking A- fe = 
(dXik/dp)p' , and taking it- along 



dw 



ihP^p'. 
dp 



(55) 



Since I s along fcj, eq. (|55|) implies ^ also along fej. 
Hence the mode is longitudinal. 

We now use ([50)) and the sound velocities of (14"6")) to 
eliminate p' from ((33J> and ([55]), Then JSJ multiplied by 

a; gives 



ipkiv'iC 2 L = -ujkiu'ipcfs, 
and (|55[) multiplied by w gives 

dX 



dw 



iPp^k 2 
dp . 



(56) 



(57) 



Since and v[ are along fc,-, eq. (|56|) implies that for 
the diffusive mode 



= ^4(4?/ C U,) = -Ui{c lS /cf L ), 



(58) 



which is independent of w. We interpret this as the lat- 
tice velocity being out of phase relative to the matter 
velocity v[ so that the fluid and lattice stresses cancel. 
Combining (151)1) and ([57| then yields 



w (cfs + c zl) = - ik 2 Pp 
= - ik 2 /3p 

Therefore 



2 dX 2 1 9A 
° ls d~p +ClL 



dX dP 
dw dp 



p dw 
dX dP 
dp dw 



dX dP 
dw dp 



d\ dP 
dp dw 



2&(A--P)-P&(A-P) 



(59) 



(60) 



For P a = 0, eq. fBDJ agrees with Ref. |27ffi and with 
Ref. For either a pure liquid or a pure solid, Dl — > 
0: 

• For a pure liquid, derivatives with respect to strain 
go to zero: dX/dw — > and dP/dw — > 0. Therefore 
Pl -> 0. 

• For a pure solid, derivatives with respect to density 
(at constant strain) go to zero: dX/dp — > and 
dP/dp -> 0. Therefore D L -> 0. 

If the system is not supersolid, and if the samples are 
not perfect, then it is consistent to interpret the obser- 
vations of Ref. in terms of this lattice diffusion mode. 

Substitution for the four derivatives in (|60[) from Ap- 
pendix [B] gives, to lowest order in P a /K. 



K 2 



V d 2 K 



K dK 



V 



2 dV 2 K+ %p v dV K 



zPv 



dK 
dV 



(61) 



where derivatives with respect to V are taken at constant 
(wik , er, N) . The form (|6T|) does not apply to case of a 
pure liquid (whereas (I60|) is general), because it does not 
permit P to have terms independent of strain. Recall 
that we have assumed that it is valid to expand K around 
P a = 0. If, in Dl, all other dependences on P a can be 
neglected, then (fBTj) implies that Dl ~ P 2 . 



V. LONGITUDINAL RESPONSE OF NORMAL 
SOLID 

Recall that (3 has units of D divided by pressure. 
As T — > we expect that, by the Arrhenius equation, 
D — > as exp[— A/fcsT], where A is a hopping en- 
ergy, because the hopping rate should yield such a de- 
pendence. Therefore, if the wavevector fc is replaced 
by d^ 1 , where d is a characteristic distance (the plate 
separation in Ref. l22t ). then the characteristic response 
time t ~ cj _1 ~ (Dk 2 )^ 1 ~ d 2 /D. Hence the view 
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that the experimental results of Ref. are due to a lat- 
tice diffusion mode leads to the conclusion that r varies 
as exp[A/fcsT]. Indeed, such a dependence is observed, 
with A ~ 30 mK. It would be useful to test for the pre- 
dicted (^-dependence. For instance, the present theory 
predicts that ch ang ing the plate separation in the pan- 
cake cell of Ref. [22| from 100 /im to 200 /im should yield 
a relaxation time approximately four times longer. 

This mode provides a means for vacancy flow to equi- 
librate vacancy concentrations. It is consistent with the 
observation of Ref. [22] that pressure decreases during an 
anneal, and when the system relaxes at constant temper- 
ature. We interpret this to mean that vacancies diffuse 
to or from the surface. 

We now turn to how a normal solid will respond to the 
two devices usually employed to generate a disturbance: 
a heater and a transducer. Since there are three longi- 
tudinal modes (thermal diffusion, lattice diffusion, and 
elastic waves), it would appear that there is need for an 
additional independent generator. Perhaps surface prop- 
erties introduce a new boundary condition that amounts 
to having an independent generator. For example, the 
material against the solid 4 He may cause the 4 He surface 
to prefer vacancies, as opposed to atoms. Thus the sur- 
face treatment may affect the behavior of both heaters 
and transducers. This argument applies to any two or- 
dinary solids, and there may be some for which this can 
be readily tested. Hence two macroscopically identical 
heaters or transducers made of different materials, or of 
the same material but with different surface treatment, 
would not show identical behavior. Since Vi — u, ~ for 
the temperature mode and the elastic modes, one way to 
characterize the response of a surface is in terms of Vi~vn. 
Thus (vi — Ui)/P' for a longitudinally moving transducer 
and (vi — iiij/T' for a heater would characterize differ- 
ences in the response to different surface conditions, and 
the extent to which they can generate the lattice diffusion 
mode. 



VI. SUMMARY AND CONCLUSIONS 

We have applied the Andreev-Lifshitz theory of super- 
solid dynamics to an ordinary solid with lattice defects 
- specifically, with vacancies in mind. At the thermo- 
dynamic level, this theory includes an internal pressure 
P, distinct from the applied pressure P a and the stress 
tensor A<j&. For the Andreev-Lifshitz theory this is nec- 
essary to permit a continuous variation from a super- 
solid to a superfluid. Under uniform static P a , we have 
Aife = (P — P a )Sik- For P a ^ 0, Maxwell relations im ply 
that P ~ Pf . These results are not conventional; Ref. l3ll 
does not include V as a distinct extensive thermodynamic 
variable, nor its thermodynamically conjugate variable 
P. In the present work many derivatives involving V are 
at fixed strain Wik, which is also unconventional, since 
normally one assumes that 6wu = —Sp/p^ Neverthe- 



less, the variables of Andreev and Lifshitz must be taken 
if vacancies are to be permitted. 

For an isotropic model, the normal modes were ob- 
tained. There are, as expected, two sets of propagat- 
ing transverse modes, with velocities as expected. There 
also are, as expected, a set of propagating longitudinal 
modes, but with velocities containing both solid-like and 
liquid- like contributions, and which depend upon P a . In 
addition there are two diffusive longitudinal modes: a 
well-known mode that dominantly involves temperature, 
and another mode involving lattice defects (i.e., vacan- 
cies). Our analysis of the physical nature of this mode 
shows that it is surprisingly complex. It involves the 
mass density p, the lattice velocity isj, and the mass-flow 
velocity Vi, with the fluid- like pressure P associated with 
p essentially canceling the solid-like stress A associated 
with Ui. 

In a separate work-22 we discuss the normal modes of 
the full Andreev and Lifshitz theory for a supersolid, 
which has nine variables. As Ref. ]l| established at T = 0, 
there are four pairs of propagating modes. Three pairs 
are essentially the elastic modes we have studied here, 
with a weak coupling to the superfluid. The fourth pair 
is basically a fourth sound mode, where the normal fluid 
is entrained by the lattice. These propagating modes, 
in the presence of a finite P a , and their generation by 
transducers and heaters, have been considered in Ref]43j. 
We also find a rather complex additional mode, not con- 
sidered in Ref. Q], which is diffusive^ Although the ad- 
ditional supersolid diffusive mode is similar to the nor- 
mal solid diffusive mode found in the present work (e.g., 
zero net stress, and distinct mass and lattice motion), its 
mode structure differs significantly. The supersolid dif- 
fusive mode is characterized by three velocities: v' ni , v' si , 
and Uj-, associated respectively with the normal mass, su- 
perfluid mass, and the lattice. For supersolid 4 He with 
P a < K, we find that vL > v' ni > u' v We also find that 
g' = p n v' n + p s v' s w 0. If He is a genuine supersolid, then 
this mode provides an alternate explanation for the expo- 
nential time-dependence of the pressure decay observed 
by Ref. [H 

We close with the following comment. Ref. 0] pre- 
dicted that supersolidity will occur because of quan- 
tum diffusion, a situation that occurs at such low tem- 
peratures that the relevant bulk diffusion processes are 
temperature-independent. Ref. [22] observe temperature- 
dependent relaxation; therefore their system is not at a 
low enough temperature to be in the quantum diffusive 
regime. Note that quantum spin tunneling is an estab- 
lished phenomenon, wherein the magnetic relaxation rate 
saturates at low enough temperatures 
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Appendix A: Andreev-Lifshitz Supersolid 
1. Thermodynamics 

Consider a general frame of reference, with non-zero 
supcrfiuid velocity v s and normal fluid velocity v n . Let 
u be the local displacement of the crystal sites relative 
to their equilibrium, and take the strain to be given by 
Wik = diUk- Then by thermodynamics the differential of 
the energy density e is given by 

de = Tds + Xikdwik + pdp + j s ■ dv s + v n ■ dg. (Al) 

Here Xik is an elastic tensor density (with the same units 
as pressure P), p is the chemical potential (with units of 
velocity squared) , j s = g — pv n (a requirement of Galilean 
relativity) , p = p n + p s (the sum of the normal and super- 
fluid densities), and g = p n v n + p s v s . By thermodynamic 
extensivity we also have 

e = -P + Ts + XikWik + PP + j s ■ v s + v n ■ g (A2) 

and the Gibbs-Duhem relation 

= — dP + sdT + WikdXik + pdp + v s ■ dj s + g- dv n . (A3) 

The system will be in equilibrium when the thermody- 
namic forces d{T, diXik, dip, diV n j, and dij si are all zero. 

2. Dynamics 

The thermodynamic variables e, s, Ui, p, v s , and g are 
taken to satisfy equations of motion that are first order in 
time and that satisfy appropriate properties under space 
rotation and inversion, and under time-reversal. Thus e, 
p, and g satisfy conservation laws (a flux but no source) , 
the phase gradient v s is proportional to a gradient (a type 
of flux, with no source), and the displacement Ui has a 
source but no flux. Thus 



d t e + diQi 
d t s + dtfi 
d t Ui 
d t v s + VO 
dtp + digi 
d t gi + d k ILik 



0, 

§, (R>o), 

Ui, 

o, 
o, 

0. 



(A4) 

(A5) 
(A6) 
(A7) 
(A8) 
(A9) 



(The source Ui was implicit in previous theoriesi^) The 
unknown fluxes Qi, fi, (f>, and Hit, and the unknown 
sources R and Ui, are determined by subjecting them to 
the condition that, when applied to the thermodynamic 
equation (|Aip . the density R of the rate of dissipated 
energy be non-negative. Note that g t is already known, 
and Qi and R will not be needed. For / i; Ui, 9, and 11^ 
we have, when terms non-linear in velocities and strains 



are neglected, 

fi SV n i 



^-djT — ^diXi k , 



= >'„, +^-d j T + p ij d l X lk , 



9 = p - Cikd k v ni - xdkjsk, 

fli/c = {PSik — Xik) — Viklmd m V n l — Qkdljsl- 



(A10) 

(AH) 
(A12) 
(A13) 



In each of these equations, the last two terms are dissi- 
pative and the preceding terms are reactive. 

Refs. [3 and [23| obtain a nonlinear term in Ui , which 
may be obtained by letting Vi — > Vi — VjdjUi. On the 
other hand, Ref. [27] obtains two nonlinear terms, which 
may be obtained by letting Vi — > Vi — VjdjUi — UidjVj.— 



Appendix B: Relevant Thermodynamic Derivatives 

In what follows, the quantities {dP/dwji) p ^, 
(dP/dp) Wikja , (dX ik /dWji)p t(7 , and (dX ik /dp) Wik:tT 
are obtained in terms of P a and the elastic constants. 

(1) With to^ = u4° } = ^33, eq. © gives 



dP 



dwik 



V,s,p 



-3K d ik w{{ = ^—dik, 
aw 



(Bl) 



where dP/dw is defined for later convenience. Substitu- 
tion for from (|3"3"|) gives, to second order in P a /K, 



dP ^ K* 
dw a K 

(2) From ([29]) we have 



pi jy-*2 

2K 3 ' 



dP 
dp 



9 (0) 2 dK* 
2«n 



By 

dK* 



dp 



dp 



± K-V™ 

dp \ dv 



(B2) 



(B3) 



<7,Wik,N , 



V 2 d 2 K 



p dV 2 



Thus (IB3I) can be written as 



dP 
~dp~ 



9V 2 (0) 2 d 2 K 
~2p~ Wl1 dV^ 



(B4) 



(B5) 



To second order in P a /K, eqs. (|33]l . (|B3|l and fB5j give 



dP 
~d~P 



1 P 2 dK* 



2K 2 dp 



V 2 P 2 d 2 K 



2pK 2 dV 2 



(B6) 
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(3) From ([23]) we have 



With 



OX 



i k 



= ( K - -p v ) S ik Sji + fi v (SijSki + SkjSu) 



We now define 



OX 4 



(B7) 



(B8) 



so that 
dX ik 



dw 



dX 

— -Q^SikSji + pv {dijSki + SkjSu — 2SikSji) . 



(B9) 



(4) From ([22]) we also have 



dXik 


j dK 




2 dpv 


dp 






3 dp 




dp v 

°P w ih ,c 


H 0) + 



x (o) 

O t kW , 



(BIO) 



dK _ K*-K 

dP w tk ,a P dVw lk ^,N p 



(Bll) 



and with a similar relation for /^y, eq. (IBlOp gives 
dX lk 



dp 



K - K 2 n v - fj, v \ (o) 

P 3 P 

Py - PV ( (0) (0) 



With (US} and = = w { 3 ° 3 \ 



dX r 



dp 



K* - K 



X (°) Q 



K* -K 



dp 



OikWli 



(B13) 



where dX/dp is defined for later convenience. To second 
order in P a /K, eq. ([3"3"]) gives 



<9A 
9p 



if* 
~K 



Pa PlK* 



2pK 2 



(B14) 



